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ABSTRACT 

A sensitivity theorem by R, M. Cohn states that for linear resistive 
circuits the ratio of a fractional change in the d.c. input resistance 
to a fractional change in an internal resistance is equal to the square 
of the ratio of the current through that internal resistance to the d.c. 
input current. The theorem is extended to show the sensitivity of 
input impedance to changes in internal impedances for an arbitrary network 
at all frequencies. Equations are developed which show the relation 
between sensitivities and instantaneous power in the frequency domain. 

An extension to the time domain makes digital computer solutions 
poss ible . 
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I. INTRODUCTION 



In 1950, R. M. Cohn [Ref. 1] presented a linear resistive network 
sensitivity theorem which was new and significant in its approach to 
sensitivity. Although the theorem is restrictive in its original form, 
it is, in fact, a particular case of a more general result involving 
circuit sensitivities as shall be discussed here. 

In Chapter 2, a proof of Cohn*s Theorem is presented. In Chapter 3, 
the theorem is extended to an arbitrary input impedance. In Chapter 4, 
the results of the previous chapter, which were developed in the frequency 
domain, are interpreted in terms of time domain responses. A general 
relationship between sensitivity and power relationships is developed. 

The results are then illustrated by means of an example for which a 
computer program has been written to obtain the appropriate waveforms. 
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II. PROOF OF COHN'S THEOREM 



Consider a two terminal resistive network consisting of branches 



hi jbs b^^, having independendent variable resistors, , ^3 , . . . . , r^ . 

The input resistance R (see Figure 1) is then a rational function of the 

r's. Consider also, that a unit curre*nt is applied to the input terminals. 

It can then be shown that R has the following property: 

P 



— = • ^ 

3r ^m 
m 



and 



5R/R 

3r /r 
m m 



1 r 
m m 

R 



m 



( 1 ) 



R 



where i is the current flowing in b . Equation (1), which is Cohn’s 
m m 

Theorem [R.ef . l], states that the ratio of the fractional change in the 

input resistance to the fractional change in an internal resistance is 

proportional to the ratio of the power dissipated in the internal resistance 

to the input power. The proof that follows is based on that of Cohn. 

Let ij denote the current flowing in, and v^ the voltage across, branch 

b.. Let the network consist of q nodes designated by n 3 ^,n 2 ,...,n with the 

J H 

current source connected from n^ to n^ . If the loops are designated by 
li,l 3 ,...,l^ and the direction of current flow, in each branch and in each 
loop assigned arbitrarily, then Kirchoff’s equations may be stated as: 



t 

; b. .V . = 0 

^ ij j 

j=i 



i = 1 , 2 , ,t 



( 2 ) 



/ 3 . . i . = g . 
tj J J 

j = l 



i = 1 , 2 , ,q 



(3) 



b. . 

ij 



a . . 



1 if the branch current and loop current have the same direction, 

“1 if the branch current and loop current have opposite directions, 

0 if the branch is not in the loop. 

1 if the branch current is directed away from the node, 

-1 if the branch current is directed toward the node. 
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= 0 if the branch is not connected to the node. 



= 1 for j = 1 
= -1 for j = q 
= 0 elsewhere. 

Define e. = i.r to be the voltage across b . in a positive direction 
J J J J 

and let the voltage of node n. with respect to node n. be E. . = f.e. with 

1 J J J 

f. = 1 if i, goes from n. to n. and f. = -1 if i goes from n to n.. If 
J J ^ J J j i ^ 



node q is grounded, = E^ = R. 



(4) 



Since Kirchoff*s laws imply the conservation of power, the povjer input 
must equal the sum of the powers consumed by the individual branches. 

This may be written as: 



E I = y i^r . = R 
1 ^ J J 

j=l 



(5) 



If equation (3) is differentiated with respect to r , it is found 

m 



that 



Br 



m 



n n 

I ) a . .1 ."I = > a. . 
L ^ ij jj ij 

j=l j=l 



Bi . 

J 

Br 

m 




Br 

m 



0 



i = 1,2, ... ,q (6) 



Let the nodes of b . be denoted by n and n , such that e is positive 

J ja jb . • j ‘ 



going from n. 



to n ,, 
jb 



then 



.2 

1 .r . = e .1 . 
J J J J 






(7) 



n 



/ 



.2 

1 .r 

J j 



j=l 



n 



- I <"ib 

j=l 





( 8 ) 



By defining a function h. (an arbitrary function of the r's) such that 



n 

a h . = 0 
ij J 

j=l 



= 1 , 2 ,. 



(9) 
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and using equation (7) above, 



n 

V" 



n 



) i.r.h. = / e .h , = /(E., - E. )h . 
J J J ^ J J Jb ja j 

j = l j=l j = l 



( 10 ) 



Since E^^ and E^^ are voltages measured with respect to node q, they 

may be written as E^'s. Also, a network is not necessarily a series 

connection of resistors. Thus, E and E. may be associated with more 

Jb ja 

than one branch. If the sum on the right hand side of equation (8) is 
performed and all terms with the same E^ are combined, it is seen that: 



n 



q n 



; i.r.h. = ; (E - E ) h. = ; ( )a. .h. )E. 

^ J J J jb ja j Z_, \ Z_. ij j/ 1 

j=l j=l i=l j=l 

'T' 

= 0 by the definition of / a. .h . 

^ ij J 

This can easily be shown in a simple example (see Figure 2). 

i - >=3>"2 + 'h - '=2>'’5 ^ 

k=l 

(E 2 - (’'S - ^>"^3 

= E,(h, + h., + hj + E^Ch^ - h^') + E.,(-h., + h., - h,) + 



( 11 ) 



1 ' 1 
E^(-h^ - ^3) 



( 12 ) 



From the definition for a. . given earlier: 

ij 

^11 "" ^12 " ^15 " ^24 " ^33 " 

®25 " ^34 " ^32 " ^ql " ^q3 " 
and all others equal 0, 

If each h in equation (12) is multiplied by its appropriate a. . 

j 

such that i corresponds to the appropriate index of the E^ , it vjill be 
seen that equation (11) holds. 



Comparing equation (6) and the defining equation for tu , (9), it is 

seen that 3i./9r = h . Thus from (11) 

J m j 



8 



0 



(13) 



n n - . 

) i ,r .h . = ) r . i . — = 

^ J J J J J &r 

j=l j=l 

If equation (5) is now differentiated with respect to r , the result is: 

m 

n n ^ n 



Br 



m 



m 



j = l 



From equation (13) 



3i 



n 

V' 2 


Sr ^ 


= ; i 


in’ -+2)1 


j 


Sr t—i 


j=i 


” j=l 


0 





m 



2 y I.r. — j 
j = l 

Since the resistors of the circuit are independent, 
Sr j- 0, j 7^ m 

Sr 1 I, j = 



m 



m 



and equation (14) becomes 



(14) 



M = -2 

dr m 

m 



(14a) 



This is Cohn's theorem for a resistive network with a unit step input 
current . 
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III. EXTCNSION_OF COH^^ * S TI-IEOREM TO AN ARBITR^\RY IiMPEDAN CE 



The expression developed in equation (lAa) is interesting, but is 
limited to purely resistive circuits. Therefore, an extension of equation 
(14a) to an arbitrary impedance is desireable. By applying a theorem of 
Tellegen [Ref. 2], it can be shown that: 



^Z(s) ^ 

9Z^(s) 



-.I(s) - 



(15) 



The method used follows to some extent the work of Penfield [Ref. 3]. 

Tellegen’s theorem states, *‘In a network configuration, imagine 

branch currents i such that for every node Si = 0, imagine branch voltages 

V such that for every mesh Sv = 0, and for every branch let the positive 

direction of current be from the + to the - denoting the positive polarity 

of the voltage. Then Siv = 0, where the summation is over all branches”. 

Tellegen then goes on to show how for Figure (1) the theorem can be 

written as SI v = S i, v, where i and v are the terminal voltages and 
Ltti^bb t t 

currents and i and v are branch currents and voltages. It is significant 
b b 

to note that Tellegen *s Theorem is strictly topological and is related in 

no way to the linearity or nonlinearity of the circuit components. For a 

constant, linear system, Tellegen shows that Si, V/ = SI,'V; . The unprimed 

b b • b b 

quantities are the first state values and the primed quantities are the 
second state values of two circuits with the same topology.^ 

Since the Laplace operator is a linear operator, it can be applied 
and Tellegen ‘s theorem can be written as: 



^It is unfortunate that Penfield chose to use the term ‘'state values” 
since it can easily be confused with the "state variables” often used in 
circuit analysis. As used here, "state values” refer to the values of the 
voltages and currents of the branches of two different circuits. These two 
circuits have the same topology but not necessarily the same circuit 
elements . 
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(16) 



or 



Vi (s)v'(s) 

L C 

t 



= ^Ib(s)V-(s) 



Vl-(s)Vb(s) =7Vb(s)Ij;(s) 

t b 



Subtracting (17) from (16) yields: 

V- r nr 






(17) 



(18) 



Unless otherwise stated, all calculation in the remainder of this 
paper will be in the Laplace domain (s - domain) . 

Consider circuits as shown in Figure 3. 



For circuit (a) For circuit (b) 



V(s) I(s) 


V(s) 


I(s) 


Port E (s) I(s) 

C 


Eg(s) 


I(s) + 61(a) 


o 

Circuit V, (s) I, (s) 

b b 


Vb(s) +6Vb(s) 


I (s) + 61 (s) 
6 b 


Using equation (18) , 


I(s)E (s) - [I(s) +6I(s)]E (s) =r 

o o ^ 


(lb(s)[Vb(s) + 6Vb 


(s)] - 



Vb(s)[lb(s) + 6Ib(s)]} 

Eg(s)I(s) = X 
b 

V- (s) = Z (s)I, (s) and for small changes 
b b b 

6Vb(s) = Zb(s)6Ib(s) + Ib(s)6Zb(s) . • (20) 

Multiplying (20) by Ij^(s) and rearranging will yield: 

I (s)6\' (s) - V,(s)6I (s) = F(s) 5Z (s) . (21) 

b b b b b b 

From the two states in Figure 3, it can be seen that: 

E (s) = I(s)Z(s) (22) 

8 

and E (s) = [l(s) + 6l(s)][Z(s) + 6Z(s)] 

6 

= I(s)Z(s) + I(s)6Z(s) + Z(s)6I(s) + 61(s)6Z(s) (23) 
Substituting equation (22) into (23) shows that: 

I(s)6Z(s) + Z(s)5l(s) + 6l(s)6Z(s) = 0 (24) 
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I (s)6Z(s) = -E (s)6I(s) - I(s)6I(s)6Z(s) 

s 

= 6l(s)[E (s) + I(s)6Z(s)] . . 

S 

If I(s)6Z(s) « E (s) [i.e. 6Z(s) « Z(s)], then; 

8 

I^(s)6Z(s) = -E (s)6I(s) 
g 

Substituting (26) and (21) into (19), it is found that: 

I^(s)6Z(s) = ^ I^(s)6Zj^(s) . 

b 

Since Z(s) is a function of Z (s) , 6Z(s) = / 6Z (s) 

b Lj is; b 

b 

Putting this result into (27) results in 



5Z (s) . ^ 

b b 

.. rl. (■s)n2 

dZ (s) " L Ll(s) J 

b ^ b 



52. 



5z(s) ^ 

azj^(s) Li(s) . 



( 25 ) 



(26) 



(27) 



(28) 



This is Cohn's theorem in the s - domain for any impedance. 

If s is set equal to jco, the results can be obtained in the frequency 
doma in : 



Z^(jcu) Ll(ju)) J 



(29) 



The right hand side of (29) may be written as a transfer function as a 
magnitude and a phase angle. 



I(j^) 

,(j^) Li(j'^) 






(30) 



(31) 



If Z|^(ju)) = R, 



_ t2 



3R 



T (gj) /29(u)) 



(32) 
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( 

. 



1 , 



If = jtuL, 

loic) - 

and 



5Z(.ia)) 

BL 



= ju)T (o)) /29 (o)) = ojT (w) /29 (ai) + r/2 



If Z^(jco) 



and 



= 1/jioC. 



iTujS) = 



11^ = (l/iuo)T^(w) /29(gj) = /29(u)) - n/2 



0) 



Thus, the magnitude of the sensitivity changes as the product of 

2 

or T (o)) /o) and its phase is advanced or retarded by tt/2 when the 
of interest is an inductor or a capacitor rather than a resistor 



(33) 



. (34) 

yjT^(oi) 

element 
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IV . COHN'S THEOREM IN THE TIME DOMAIN 



Equation (28) expresses Cohn's theorem in the s - domain for any 



input current I(s) applied to an impedance Z(s). 



SZ 



_ r 

^(s) Ll(s) J 



This equation is easy to use if the input current I, and the internal 

element current are known as functions of s. The input current wave 

shape v;ill be known for most application and its transform to the s - domain 

will usually not be difficult. However, for any but the most elementary 

of circuits, the internal element current can be found only through 

a series of tedious calculations. In contrast, there are several good 

computer programs which will calculate i (t) . It is, therefore, desireable 

b 

to find a direct means of determining the sensitivity as a function of 
time rather than by taking the inverse transform of (28). For the class 
of input currents given by (39), as derived later, it can be shown that: 

2m 



SZ(s) 

BZ (s) 
b 




k-2n^2 



(s) 



k=0 



(35) 



where n and m 
results in the 



are integers 



time domain 



2n 

5 z __ Y 

dx ' 

b 



defined by 



^k-2n+j 
k ^^k~2n+j 



(39). This leads 




(t) .V i^(t)] 



to the following 



(36) 



where the symbol * is used to indicate convolution, and 
j = +1 if is an inductor (L) 

= 0 if x^ is a resistor (R) 

= “1 ifx^is 1/capacitor ( I /C) 

A very interesting simplification occurs when i (t) = 6(t), a unit impulse 
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i 



so that E(s) =Z(s). 




= \(s)Ib(s) 



(37) 



(38) 



Equation (37) shows that the sensitivity is related to the instantaneous 
power in the frequency domain. 

In attempting to develop the time domain equivalent of equation (35), 
it would appear to a casual observer that a simple application of the 
principal of convolution would apply. That is, X(s)Y(s) = ( t) v-y ( t ) ] 

and thus, X(s)/Y(s) = •^[x(t) vrVy(t)]. However, this is not true since 



be convolved with x(t) . There does not appear to be any such simple 
operation . 

An alternate possibility is to limit the consideration to input 
current wave shapes which can be manipulated easily so as to yield a direct 
solution. Consider an input current of the following form: 



^(t)] 7^ Y ^(s) . For example: 



Let Y(s) = s 



y(t) = a unit doublet 



Now let 1/Y(s) = 1/s 



1/s transforms to a unit step in the 



time domain and a unit step is 



not the reciprocal of the unit 



doub let . 



It would appear useful to try to determine if there is a direct way 



to obtain the inverse transform of 1/Y(s) to get a time function which can 




( 39 ) 



15 



g = any real number 



n = any integer , or 0 

= all real numbers (all but one can be zero). 

This allovjs the input to include impulses, steps, ramps, sinewaves, decaying 

. , sin at sinh at 

or increasing exponentials, , , etc. 

a a 

Considering the input current specified in equation (39), equation (28) 



becomes 



az(s> ” 



SZ^(s) 

b 



N^(s) 



(40) 



2m 

-2 -2n ^ k 2 

= 8 s 



k=0 



2m 



V -2 k-2nV2 , ^ 



k=0 

2m 



= I( 






k-2n\^2. . 



k=0 



(41) 



The result is equation (35). 

Consider the impedance of the elements. If the element is an 
inductor, then = sL and equation (28) may be written as 







a(sLj^) ( 


VI (s) ) 


SZ(s) _ / 




* 


^I(s) / 



(42) 



If the element is a capacitor, then Z = l-/sC and equation (28) becomes 

b b 



- S Z(.s ) ._ _ 1 
B(l/C^) s \l(s) J 



(43) 



^The subscripts are changed because the b*s are now the coefficients 
of the squared polynomial. 
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If the element is a resistor, equation (28) is 

5_Z_Cs.). ^ ^VlV 
dR^ V I(s)/ 



( 44 ) 



Thus in general the sensitivity can be written as 

ifVl 

Vl(s) 



SZ(s) _ 

BXb 



(45) 



where i = +1 if x, is an inductor 
b 

= 0 if X, is a resistor 
b 

= -1 if X, is a compliance 
b 

Combining (45) and (41) yields the result 

2m 



= j 
s"b 



2m 



i( 

k=0 






k-2„)j2^^) 






(46) 



k=0 



It can be shown that multiplication by s in the s - domain is equi- 
valent to differentiation with respect to t in the time domain, that 
division by s in the s - domain is equiva lent • to integration in the time 
domain, and that multiplication of F(s) in the s - domain is equivalent to 
convolution in the time domain. Thus, if the exponent of s in a term of 
equation (41) is greater than zero, the individual terms of (46) become 

k-2n p 

i, (t) * i (t) ! 



k k“2n L. b 
o t 



J 



If the exponent of s is negative, the term of (46) becomes 

* 1 ^( 0 ] . 



'Sec, for instance, Holbrook [Ref. 4] pp. 66, 68 and 108. 
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The number of integrations is equal to k-2n. The sensitivity, 



CbZ 



then becomes a sum of integrations and/or differentiations of the 
convolution of the current through the element of interest uith itself, 
and equation (46) can be v?ritten as: 

2m 



Bz(t). ^ Y , 

sx, L 

^ k=0 



k-2n+j p 
k ^^k-2n+j 



This is equation (36) . 

To show the power relations of equation (37), Ohm's lav? is applied 
(see Figure 4). Again, unless otherwise stated, all computations are made 
in the s - domain. 



5Z(s) 

SZ^(s) 






Thus : 



E (s ) = I ( s)Z(s) 
g 



SE(s) = I 3Z(s) + Z(s)dl (s) 
g g 



SZ(s) 



= 5E(s) 

I (s) 
g 



5E(s)/I (s) 

S 

(s) 

b 



rl, (s)n 

L^J 



2 



SE(s) ^ ^b^^^ 

BZj^(s) Ig(s) 



BZb( 

If the input current is limited to an impulse in the time domain, I (s) 

S 

will be unity and equation (42) becomes 



^H(s) 

s)/Z^(s) 



I,(S) 



I,(S) 



Ig(s) 



Eb(s)Ib(s) 

Ig(s) 



(48) 



5ECsj., „ _ 

5Z^(s)/Z^(s) 



= Ej^(s)I^(s) 



(49) 
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If a circuit is excited with an impulse current (see Figure 5) , the 
ratio of the voltage across any element to the input current may be 
written as a transfer function such that: 



I^(s) 



= T(s) 



(50) 



E(s) = T(s)l^(s) 



If the circuit is now excited by a current source equal to the current 
through the element of interest, the new element voltage, E(s) , will be 

(51) 

This is true since the transfer function T(s) is a function of the network 
and not of the input current or voltage. The original input current was 
an impulse so that equation (50) can be written as E^(s) = T(s) . 



Substituting this into (51) will yield 
E(s) = Ej^(s)I^(s) 



(52) 



The input current can be vjritten in the time domain as i. (t) , and E, (s)I. (s) 

b b b 

can be written as e^(t) * i^(t). Thus, (52) becomes in the time domain 



c(t) = ej^(t) * 



(53) 



Equation (53) is the right hand side of equation (38). Thus, in order to 
measure the sensitivity of equation (38), the following procedure can be 
used. Excite the original circuit by an impulse and calculate i^(t). 

Then excite the circuit with i^(t). The voltage across the element is then 
equal to the sensitivity. 

There are several interesting variations of equation (28). For 
instance: since Y(s) = Z ^(s) 



3Y(s) ^ 1 5Z(s) ^ 

■ z2(s) ■ ■; 

8V(s) . _ 

BZ^(s) L E(s)j 



-I (s)^2 



T (s) -,2 



Li(s) J Li(s)z(s)J 



(55) 
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I 



(56) 






SZ^(s)/Z^(s) “ '^b^^^LE(s) J 



,-Ij^(s)^2 ,-E, (s)^2 



b' ' b 

If I = EY is substituted into equation (28) 



, , -T 

LE(s) j 



5Z(s) ^ 

BZ^(s) 



I (s)_2 _E^(s)Y^(s)_,2 



Ll(s) 



.1 =r_b:_:jLiii 

J Le(s) Y(s) j 



v^s)5Z(s) ^ 

Yf(s)SZ (s) 
b b 

- 1 2 
Remembering that Z = Y , the partial derivative yields SZ =-5Y/Y 

Putting this into (56) results in 

5Y(s) ^ 

3Y^(s) L E(s). 

For another variation, Z = Y ^ is used. Thus 



(57) 



5Z(s) ^ 

SYb(s) 



Y (s) 

5Z(s) ^ 

5Yj^(s) L I(s)J 

52 (s) . 

5Y, (s)Y (s) 
b b 



1 5Y(s) 1 

.2. . 5Y, (s) ■ 2, , Le(s) j 



Y (s) 



Yb(s)Eb(s) 

I^(s) 



I^s) 



(58) 

(59) 



If e(t) is restricted to a delta function, then equation (55), (56), and 
(57) may be written as 



5Y(s) 

SZj^(s) 



-I^(s) 



5Y(s) 

5Zj,(s)/Z^(s) 



-Eb(s)Ij^(s) 



(60) 



(61) 



5Y(s) 

SY^(s) 



= 



(62) 



If i(t) is restricted to a delta, equation (58) and (59) may be written as 
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i 








( 63 ) 




5Z(s) 



= -I,(s)E^(s) 



(64) 



SYj^(s)/Yj^(s) 



If the denominator of the right hand side of these equations are 
restricted to the forms of equation (39), time domain functions similar to 
that of equation (47) can be obtained. For any sensitivity calculation 
however, the impulse driving function gives the simplest and most direct 
result . 
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T\vo Terminal Netv;ork vith Variable r, 

k 

Used for Development of Cohn's Theorem 
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Figure 2 



Five Branch Network Used to 
Demonstrate Cohn's Theorem. 
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a) State 1 




Z -f 6Z 

b) State 2 



Figure 3 

Circuits with the Same Topology 
But Vv7i th Different Internal Impedances. 



2A 




Figure 4 

T\vO Terminal Networks Used for Developing 
the Power Relations. 

If a constant input current is required, a 
current source is connected to the input 
terminals. If a constant voltage source is 
required, a voltage source is connected. 

In either case E = IZ . 
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i 



Ig(s) 



T(s) 





Figure 5 

Two Terminal Network Used for Developing 
the Time Domain Solution. 
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APPENDIX A 



Example Using a Parallel RLC Circuit 

-4 -7 

For the circuit in Figure 6; let R = 1.0, L = 10 , C = 10 

Calculate the sensitivity of the input impedance with respect to each of the 
circuit components for a delta input current. Also, determine the sensitivity 
of the input voltage with respect to 6Z/Z for each component with the same 
delta input. 

An impulse was approximated by a triangular pulse 8.0(10 seconds wide 
and 2.5(10^) amps at its center. This input current was used as data for 
an ECAP solution to arrive at the currents through and voltages across the 
components . 

The data obtained from the ECAP solution was used as input data to 
Subroutine SEN which is in the Computer Program section of this paper. Since 
outputs were obtained from ECAP every 50 calculations, the time increment 
for SEN was 5(10 ^) . The final time was 5(10 ^) . The delta was chosen 
to give it unit area so COEF is 1.0. The Laplace transform of a delta is 
a constant. Thus, SN and IPR were set at zero. The lELEM and ITYPE* values 
were determined by the type of element and type of sensitivity required. 

The results are shown in Figures? through 12. 
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Figure 6 

Parallel RLC Circuit Used for Sample Calculations. 
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Figure 7 - Change in the input impedance with respect 

to a change in an internal inductance. 
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Figure 8 - Change in the input impedance with respect to a change in an internal 
resistance . 
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Figure 9 - Change in the input impedance with respect to 
change, in an internal reciprocal capacitance. 
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Figure 10 - Change in the input voltage with respect to a fractional 

change in reciprocal capacitance for an impulse input current. 
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Figure 11 - Change in the input voltage with respect to a fractional 
change in resistance for an impulse input current. 
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Figure 12 - Change in the input voltage with respect to a fractional change 
in inductance for an impulse input current. 
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PROGRAM 



CO'^PUTPR 



SUBROUTINE SON ( H , T F , C JE P f S N , I PR , K E Y 1 , 1 c L cM , I T Y P E , A-, 
1B,F,T, SFIIS) 

PURPOSE 

TO COMPUTE THE SENSITIVITY OF THE 10^1)'’' IMPCOAMCF 
OR INPUT AOMITTiNCE 0^ AMV RLC CIPCHIT WITH RESPECT TO 
1 , THE CHANG:: IN ANY IMT:"-;i\AL J M'->roA'iC T 

P. TI^l CHANG.: IN ANY INT''RMAL A I) ''O Y T /\NC F 

i. Tit:: chang: in /hny i mp.. ):\n('-' .oiviivji) by that 

1 MPIMaNC.P 

.'rc Tl-F ^'HANr.. IN an AU'UTTAN^’’ 'OlVUn'O BY THAT 

AO MITTANC: 

THE SUBROUTINE PIIL ALS1 CALUJLATE T(i): S’NN S I T I V I ty OF 
TI'P INPUT VOLTAG- wRT iX-L'^A / nvr- R 7 . , wupR,' / is ANY 
lMTi'Rf''''.L I MPii-JC \CH . the -L -CTR IC.O. SOOKC;-; ^UlST B i 
LXPIO.SSEI.) AS A LA PL’ ACE TIv A N S •- "i;? M . 



USAGE 

CALL S!.\( IT ,TF ,C 

T.if) oat?. Cards 
GRAPH i'JTPUT IR 
THE FJHSr “T C'H 
CARO TH.; SECOND 



■'vF ,SM, I P-' , K: Y1 , 1 Tl L'^ , 1 TYPE, A, B, E, SENS ) 

ARE RvOUI ’.’.M F0;> TIN-: GRAPH TITLE 1^ A 
RFP:Ui-ST,-i^. TH- FIRST CARD It 1ST CHMTAIN 

;pAC t: p ■> ^ ihl titL(. a'to the sfooino 



DESCRIPTION OF PARAMfcT'^RS 

H - THE INCR = ’^''FNT OF TH= ARGUMENT VALUI^S 

TF - THE PINAL TIV- OF TH?- INPUT VECTOR (A) 



COEF 

SN 

IPR 

KEY! 

lELEM 



ITYPE 



A 

B 

E 

T 

SENS 



VALUE OF THE SOURCF NU^'^PATOR CO-FFICIFNT 
POWER OF TH<^ S TFRM T(!E SMUPCP MU'^FRATOR 

HIGHEST POWER IN THE SOIIRCR NUNEPATPR 
CONTROLS graph OUTPUT: 0 PROPHCSS NO GRAPH 
lOEMTIFIFS THE ELPMPMT OF INT^^RFST; ENTFR A 1 

fcr an indhcto'^, a 0 for a rfsistop, or a -1 

FOR A CAPACITOR 

TYPE OF SENSITIVITY REQUIRRD. ENTER A 0 IF 
THE CONVOLUTION IS TO BE BETWFFN DIFFtpcht 
VECTORS. ENTFR A 1 IF THF CONVOLUTION IS TO B 
BETWEEN IDENTICAL VFCTHPS. 

CURRENT OR VOLTAGF VECTOR FOR th= ELE'-'ENT QP 
INTEREST 

CURRENT OF VOLTAGE VECTOR FOP THE ELEMENT PF 
INTEREST. 

COEFFICIENTS OF THF SOURCE DENOMIKiatOP IN 
ORDER OF I NCR easing P0W“R. 

OUTPUT VECTOR 0^ THE T I ■''f 
OUTPUT VECTOR OF THF SENSITIVITY 



REMARKS 

1. THE EQUIDISTANT VECTOR OF THE ELEMENT CHRRfmj or 
VOLTAGE(A) MUST BF FijRNISHFO WHEN THE S^riSTTIVITY 
SOLUTION REQUIRES THE SELF CONVOLUTION OF CURRENT 
OP voltage, the EOHI distant VECTOR (B) ‘OiS-r BE 
FURNISmfD also but may PE ZERO UNLESS th^ 
SEMSTTIVITY SOLUTION REQUIRES THf CONVOLUTION OF 
DIFFERENT VECTORS. BOTH VFCTQRS MUST BE IN B515.7 
FORMAT. 

2. MAXIMUM length QF (A) AND (B) IS 100 POINTS. 



SUBROUTINES AMD FUNCTION SuppRQQRAMS RFOUIRFD 
DAHC02 {H,TF1 ,tf2,A,B,D) 

DET3 ( H,Y,7,N0I M, I<=R) 

QSF (H,Y,Z,NDIm) 



35 



ooo r>oo ooo ooo ooo o ooooooooooooo 



METHOD 

SQUARES A P'^LYNDMIAL, EXTRACTS THE POUF.RS AND US'^S 
THESE POWERS TO DET'=RM|N!= J\i'‘ NUMRER pc p)! FF £P FMT I A- 
TIONS OR INTEGRATIONS TO BE PEPccRMFD. 

REFERENCES 

(1) JAMFS G. HOLBROOK, LAPLACE TRANSFORMS Fpo 
ELECTRONIC ENGINEERS, PrRr,Ai.*ON dpcjs, 19AA. 

(2) A.L. DAHLVIGtAPPLICATinNS ^F COmm'S SENSITIVITY • 
THEOREM TO TTf-'= DOMAIN RESPONSES, ^ASTERS THESIS, 
NAVAL POSTGRADUATE SCHOOL, 1970. 

DP’^ENSIGN A ( 100) ,B( 100) ,C( 100,100) , D ( 201 , 2 ) , Y ( 201 ) , 

IZI 201 ) , W( 201 ) , SENS! 201 ) ,E( 11 ) ,0 = N( 25 ) , T(20l ) 

INTEGER SD(2'5),SN 
REAL LARl/^H / 

REAL’i'B ITITL(12) 

I2=TE/h 
13 = 12+1 2+1 
SN=SN+SN 
IP1=IPR+1 
L=IPR+IPR+1 

determine the type of SENSITIVITY REQUIRED 

IF ( ITYPE.EO.O) GO TO 81 
DO 84 1=1,12 
B( I )=A( I ) 

84 CONTINUE 
81 CONTINUE 

INITIALIZE DEN AND SD 

SD( 3 )=0 
0£N( 1)=C.000 
DO 11 I=2,L 

OEN( I ) =0.000 
SD ( 1 ) = I - 1 

11 CONTINUE 

SQUARE THE INPUT POLVNO'^'IAL 

DO 1.2 1 = 1, IPl 
DO 13 K=1,IP1 
J=I+K-1 
DO = E( I I’SEI K ) 

DL- NM J ) =DO + DEN( J > 

13 CONTINUE 

12 CONTINUE 
COEF=COtF«COEF 
DO 4 I=1,L 

IF(DEM( I ).E0.0.0000)GO TO 14 
GO TO 4 
I A SD(I) = 0 
4 CONTI NLF 

WRITE ( 6,906) 

WRITE! 6,111) ( I,DFN(I ) ,SD( I ) ,1=1,L) 

PREPARE FOR CONVOLUTION 
TFZ=TF 

CONVOLVE INPUT VECTORS; INI TALI 7E SENS 

CALL DAHC02 ( H , TF , TF Z , A , R , D ) 

DO 16 1=1,12 
T( I )=D( I , 1 ) 

Y( I )=D( I ,2) 

SENS! I ) =0.00000 
16 CONTINUE 
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DO 21 M=1,L 
J=SD( M)-SN 



DETERMINf} ELEi'^^SNT TYPE AMO SET VALUE OF J 

IF ( TELE‘S) 76,75,76 

74 J=J-L 

GO TO 75 

76 J = J<-1 

DETERMINE UHETH'^R INTEGRATION, DIFFERENTIATION OR 
NO ACTION IS REQUIRED 

75 IF (J) 22,23,24 

BEGIN REQUIRFD INTEGRATIONS AND EVALUATE SENS 

22 J=-1=?'J 

CALL OSF(H,y, W,I2) 

WRI TEC 5,907) 

IF (J.LT.2) GO TO 73 

DO 31 K=2,J 

CALL 0SF(H,W,W, 12) 

WRITF(6,907) 

31 CONTINUE 
GO TO 73 

73 no 41 1=1,12 

SENSC I )=SENS { ! )r (DENC M) /COFD^UC I ) 

41 CONTINUE 
GO TO 21 

NO ACTION IS REQUIRED, EVALUATE SENS 

23 DO 42 1=1,13 

SENSC I )=SENSC I I + COENCM) /C0EF)*yC I ) 

42 CONTINUE 
GO TO 21 

BEGIN REQUIRED DIFFERENTIATIONS AND EVALUATE SENS 

24 CALL DET3CH,Y,W,I2,IER) 

WRITE C6,104) lER 

I FC J.LT.2) GO TO 99 
DO 32 K=2,J 

CALL DET3CH,W,W,I2,IEP) 

WRITE C6,104) lER 

32 CONTINUE 

99 DO 43 1=1,12 

SENSC I ) = SENSC I ) + C DE N ( N ) /COE P ) >!' W C I ) 

43 CONTINUE 
21 CONTINUE 

determine element type for OUTPUT PURPOSES 

IF C IELFM)15,77,78 
15 WR I TEC 6,909) 

GO TO 64 

77 WRITE C6,910) 

GO TO 64 

78 WRITE C6,903) 

64 IF (KEYl.EQ.O) GO TO 79 

READ C 5,911) C ITTTL C I ) ,1=1,12) 

CALL DPAVv'C 12 ,T ,SENS ,0 ,0, LABI , I TTTL , 0,0 ,l , 1 , 1 , 1 , 6, 4, 1 , 
ILAST ) 

79 CONTINUE 

101 FORMAT C 5F15. 7) 

104 FORMATCPX,' lER = *,1^) 

111 FORMAT C 3X,I3,3X,-15.7,6X,I3) 

905 FORMAT C 3X, ' TERM ' ,5 X, » COEFE ICI ENT ' , 7X, ‘ POWER • ) 

907 FCRMATC 1H0,5X, 'SPLINT' ) 

90S format C IHO ,5X ,' TME EicmfNT IS AN IMHiiCTOP') 

909 FORMAT C IHO ,5X, ' THE ELEMENT IS A CAPACITOR') 
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910 format I1H0,5X» ‘THE ELFMcmt IS A RESISTOR') 

911 FORMAT (6A8) 

RETURN 

END 



SUBROUTINE 0AHC02 ( H , TF 1 , TF ? , A , B , D ) 

DIMENSION A( 100) ,B( 100) ,C( 100,100) ,0(201,2) 

N=TF1/H 

M=TF2/H 

NM=Kii-M + l 

DO 1 I=1,NM 

1 D( I ,2) =0. 

DO 5 I = 1 , N 
00 4 J=1,M 
C( I ,.J) =A( I )*B( J) 

K=H-J 

4 0(K,2)=D(K,2)■^C( I , J)=«'H 

5 CONTINUE 

0( 1,1)=0. 

DO 7 L=2,NM 
AA = L"1 . 

7 0(L,1)=AA*H 
RETURN 
END 



SUBROUTINE DFT3 ( H , Y , 7 , NOI M, IFR) 



DIMENSION Y( 1) ,7 ( 1 ) 

TEST OF DIMENSION 
IF(NDIM-3)4,1,1 

TEST OF STEPSI7E 

1 IF(H)2,5,2 

PREPARE DIFFERENTIATION LOOP 

2 HH=.5/H 
YY=Y( NDIM-2) 

B = Y( 2 M-Y( 2 ) 

B + B~Y( 3)-Y( 1 )-Y( 1 )-Y( 1 ) ) 

START DIFFEPENTI ATI CN LOOP 
DO 3 I=3,NDIM 
A=B 

R=MH*( Y( I ) -Y( 1-2 ) ) 

3 Z( 1-2 ) =A 

END OF DIFFFRENTIATION LOOP 

NORMAL EXIT 
IER = 0 

A= Y( NOI M-l ) +Y( N^IM-1 ) 

Z( NDI M ) =MH*( Y( NOI M) +-Y (NOI M) + Y( NCI M) - A-A<-YY ) 

Z(NDI'^-1 )--^B 

RETURN 

ERROR EXIT IN CASE NDIM IS L^SS THAN 3 

4 IER=-1 
RETURN 

ERROR EXIT IN CASE OF ZFRO STEPSI7S 

5 IER=1 
RETURN 
END 
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SUBROUTINE OSF ( H , Y , 7. , NDI MJ 
DIMENSION Y( 1) ,7(1) 



HT = . 33 33?3 3<^H 
IF(N0IM-E)7,8,1 

NOIM IS GREATER THAN 5. PREPARATIONS OF INTEGRATION 
LOOP 

1 SUM1=Y( 2)+Y( ?) 

SUM1 = SLM1 + SIIM1 

SUM! =HTY( Y( 1 ) +SitMH.Y( 3 ) ) 

AUX1=Y(A)+Y{4) 

AIJX1 = ALX1+AUX1 

AUXl^SU^'l +HT«( Y( 3) + AUX1 + Y{ - ) ) 

AUX2 = HT*( Y( 1 ) + 3. 875*{ Y( 2 1 + Y( = ) H-2.625*{ Y(3 )+Y(-£v) ) + 
1Y( 6 j ) 

SUM2=Y( 5 ) fY{ 5 ) 

SUM2=SLM2+SUM2 

SIJM2 = AUX2-HT*{ Y(4) +-SII^^2<-Y ( A ) ) 

7 ( 1 ) = 0 . 

AUX=Y( 3)+Y(3) 

AUX=AUX+AUX 

7 ( 2) =SU'''?-HT>t‘( Y{ 2) f AUX + Y (A ) ) 

7( 3)=SU'^1 
7( 4) =SUM2 
IF(\iOIM-6)5,5,2 

INTEGRATICN LOOP 

2 DO 4 1=7,NDTM,2 
SUM1=AUX1 
SUM2=AUX? 

AUX1=Y{ I-l)+Y{ I-l) 

AUX1=AIX1+ AUXl 

AUXl = SLMli-HT*(Y(I-2) + AUXUY( I ) ) 

7( I-2)=SUP1 
IF{ I-NDIM) 3,7', 6 

3 AUX2 = Y( I ) ^-Y( I ) 

AUX2=A LX2^-AUX? 

AUX2=SLM2 + HT>M Y( I- I ) +AUX2 + Y ( T + 1 ) ) 

4 Z(I-n=SU^-'2 

5 Z( NDI M-l ) xAUXl 
Z{N0TM»=^AUX2 
RETURN 

6 7( NDIM-1 )xSUf-'2 
Z{N0I^M=AUX1 
RETURN 

END OF INTEGRATION LOOP 

7 IF(N0IM-3)12,1 1,R 

NOIM IS FOUAL TO 4 OP 5 

8 SUM2 = 1.1?5=<=HT*( Y(1M-Y{2)+Y{2H-Y{2)+Y(3) + Y{3H-Y(3) + 
1Y{4) } 

SUM1=Y{ 2)+Y(2) 

SUMI = SUW1+SU''1. 

SUM1=HTV( Y( 1)+SUM1+y(3)) 

Z( 1 )=0. 

AUX1=Y{ 3)>Y( 3) 

AUX1=ALX1+AUX1 

Z{ 2) =SUV2-HT*( Y ( 2) +AUX1+Y( 4) ) 

IF{fJOIM-5) 10,'R,Q 

9 AUX1=Y{4H-Y(4) 

AUX1 = ALX1 <-AUXI 

Z( 5)^SUM1■^HT«( Y( 3)+AUXl4-Y( 5 ) ) 

10 7{3)=SUV1 
7{ 4)=SUM2 
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RETURN 

NDIM IS EOUU TO 3 

11 SU'^l=HT=i<( 1.2 5*Y{1)+Y{2H-Y{2)-.25*Y(3) ) 
SUM2=YI 2) <-Y( 2> 

SU^■2-SUM2 + SU'''2 
Z(3)=HT*(Y< 1 H-SUM2<-Y( 3 n 
Z( 1)=0. 

7( 2)=SUN'l 

12 RETURN 
END 
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